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Growth Equation with Conservation Law 
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A growth equation with a generalized conservation law characterized by an integral kernel is 
introduced. The equation contains the Kardar-Parisi-Zhang, Sun-Guo-Grant, and Molecular-Beam 
Epitaxy growth equations as special cases and allows for a unified investigation of growth equa- 
tions. From a dynamic renormalization-group analysis critical exponents and universality classes 
are determined for growth models with a conservation law. 
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I. INTRODUCTION 

In order to describe the dynamics and scaling of in- 
terface growth, Kardar, Parisi and Zhang proposed a 
Langevin equation now known as the KPZ equation |Q. 
This equation is constructed to describe the long-time 
long-wavelength (hydrodynamic) limit of the dynamics 
of noncquilibrium interface growth processes. The KPZ 
equation has been studied intensively by analytical and 
numerical methods and a number of theoretical results 
have been obtained (2]-[3] . Specifically, the scaling of the 
correlation function C(x,t) = x 2a g(t/x z ) characterized 
by the roughness exponent a and the dynamic exponent 
z has been established. However, only poor agreement is 
obtained between the theoretical results and experiments 
Jl,!). As a result, various modifications of the KPZ equa- 
tion have been analyzed. 

In the present paper we introduce a new growth equa- 
tion with a generalized conservation law described by 
an integral kernel. We will refer to the equation as 
the growth kernel equation (GKE). It contains the KPZ 
equation as a special case. In addition, the previously 
studied Sun-Guo-Grant (SGG) j|] and Molecular-Beam- 
Epitaxy (MBE) @] equations are also contained in our 
general equation. 

The motivation for introducing the GKE equation is 
to gain information on how conservation laws change the 
universality classes for nonequilibrium growth models, 
and to allow for a unified description of growth models 
studied so far. Furthermore, one can speculate whether 
some growth experiments, which yield exponents that do 
not agree with the KPZ exponents, may contain nonlocal 
growth effects such as, e.g., the experiments on electro- 
chemical deposition reported in Refs. 

Previous studies of nonlocal terms in interface related 
topics include fluctuating lines in quenched random envi- 
ronments. Domain walls subject to quenched long-range 
correlated impurities were studied in jl0| . More recently, 
the dynamic relaxation of drifting polymersj ll[ and the 
critical dynamics of contact line depinning |L2 | were stud- 
ied, where in both cases the equation of motion includes 
nonlocal interaction terms. 



II. GROWTH KERNEL EQUATION 

The GKE equation for a d dimensional interface h(x,t) 
reads 

^ = J d d x> K(x - x 1 ) (vV l2 h + ^(V'/i) 2 ) + 77, (1) 

with an additive noise rj(x,t) whose correlations will be 
specified below. The kernel K(r) describes nonlocal in- 
teractions in the system |l3). We want the total height 
H(t) = J d d xh(x, t) to be conserved and impose the con- 
straint J d d x' K(x — x 1 ) = 0, which leads to dH/dt = 0, 
provided the noise is chosen to satisfy rj(k = 0, t) = 
(see Eq. (||) below). Consequently, the GKE equation 
conserves the quantity H(t). 
The kernel has the behavior 



K(r) 



for 



(2) 



where we have introduced an exponent a describing the 
long-distance decay. By Fourier transforming, a = cor- 
responds to the kernel being a Dirac delta function, and 
therefore the usual KPZ equation (!]]. 

The case K(x-x') — — V 2 , S d (x — x') yields the dynam- 
ics of the SGG and MBE equations HQ , and corresponds 
to a = 2. In order to incorporate these equations, we in- 
troduce a kernel N(r) in the noise correlator 



(r)(x, t)r](x',t')) = 2D N(x - x')5(t - t') 



with the form 



N{r) 



for 



(3) 



(4) 



Here, r is an exponent independent of a. The case 
t = means no correlations in the noise (KPZ, MBE) 
whereas r = 2 corresponds to conserved noise as it ap- 
pears in, e.g., the SGG equation. The correlator (||) im- 
plies that 77(A) = 0, t) = 0, as required above. We can 
imagine the noise related to a white noise £(x, t) through 
r](x,t) = / d d x' R(x - x')£(x',t), where the kernel R(r) 
has the large argument behavior R(r) ~ l/r d+T//2 . This 
form leads to the correlations in Eqs. (||) and (||). 



f 



Under the rescaling x — » x' — x/b the parameters in 
the GKE equation change as 



v 
X 

D -> D' 



A' = b z+a - 2 - a \, 

LZ — 2ot~d — T 



D. 



For A = 0, the equation is made scale invariant for 



2 + a 



and 



2 + a - d - t 



(5) 
(6) 
(7) 



(8) 



If we use these values in the rescaling for A we obtain 
that A' = &( 2+cr_d ~ T )/ 2 A, so naively we expect the criti- 
cal dimension of the model to be given by d c = 2 + a — r. 
Therefore, for d > d c the A term will scale to zero, 
whereas for d < d c the A term will be relevant and the 
scaling behavior of the GKE equation will no longer be 
described by the naive exponents. Now we will carry out 
a dynamic renormalization group (RG) analysis in order 
to determine the scaling behavior of the GKE equation. 



III. RENORMALIZATION GROUP ANALYSIS 

We Fourier transform the GKE equation using the 
rules for Fourier transformation of convolutions and 
products and obtain in the hydrodynamic limit k — ► 

h(k, u>) = G (k, uj) v (k, u ) - ~ G (k, to) k° J A 

^ q-(k-q)h(q,n)h(k-q,cj-fl), (9) 

where Go(k,u>) is the (bare) propagator defined by the 
expression Go(k,u>) = l/{vk 2+cr — iu>). A is the momen- 
tum cutoff. The noise in Fourier space takes the form 

(r](k, w)ri(k', u')) = 2D k T (2Tr) d+1 5 d (k + k')5(u + u'). 

(10) 

The renormalization group consists of coarse-graining 
followed by rescaling fl4[| . Coarse-graining: Modes with 
momenta e~ e < k < 1 (A = 1) are eliminated from the 
equation of motion. Rescaling: Wavevectors are rescaled 
according to k — > k' = bk, with b — e l . The RG proce- 
dure is most efficiently carried out by the means of dia- 
grams, i.e., we represent the GKE equation (|^) as shown 
in Fig. % cf. Refs. 

After a standard but lengthy calculation one obtains 
the one-loop RG flow for the GKE equation [[17| 



dv ( X 2 D 2 + 2a - t 
-dl = V V- 2 - a + Kd — M~ 



2 -a), 



dD 
~dl 



D z~2a-d 



K d X 2 D 



(11) 
(12) 
(13) 



where A Tt<T = 1 for r > 2a, and A Ti<7 = for 
r < 2a. The geometrical factor Kd = 5*^/(2^) d , and 
S d = 2TT d / 2 /T(d/2) is the surface area of the d dimen- 
sional unit sphere. Furthermore, the RG flow for the 
case t > 2a reduces to the case r = 2a 0- Due to the 
fact that A does not renormalize to one-loop order, cf. 
Eq. (|l2|), one obtains the exponent identity 



a + z 



(14) 



The critical exponents are determined from dvjdt = 
and dD / dl — 0, i.e., by fixing v and D. It is convenient 
to introduce the coupling constant 



9 = g{t) 



K d X 2 D 
Ad v s 



(15) 



with the dimension [L] d+r 17 2 . The RG flow of g for 
fixed v and D becomes 



da 1 dX 

i = 29 xn = 29{a+z -°- 2) 



(16) 



(2 + a - t - d)g - [3(2 + 2a - r - d) - dA Tt<J ]g 2 



The fixed points (FP) for g are 

2 + a - t - d 



9a = 0, 9 = 



3(2 + 2a - t - d) - dA Tta ' 



(17) 



Physically, we have that v, D > 0. As a result, FPs 
where g* < are unphysical since they would lead to an 
imaginary value for A. 



IV. RESULTS 

First, we discuss the case r < 2a. The critical dimen- 
sion is d c = 2 + a — t. In Fig. || we show the RG flow 
of the coupling constant g, for various dimensions. For 
d < d c the trivial FP <?q = is unstable, whereas g* is 
stable. For d > d c , g$ becomes stable. For d > d c + a, g* 
is again positive but now an unstable FP. Probably this 
latter behavior is an artifact due to the fact that the FP 
is only known from the one- loop expansion Jl7j| . 

The FPs and exponents can be calculated in an 
e = d c — d expansion. For the trivial FP = we 
recover the exponents for the linear GKE equation, cf. 
Eq. (||). Consequently, they will describe the GKE sys- 
tem for dimensions d > d c , where the trivial FP is stable, 
cf. Fig. ||(b). The nontrivial FP is to first order in e equal 
to g* = e/3cr, with the exponents 



and 



e 2 + a - t - d 

a = — = , 

3 3 



e d + r + 2a + 4 
Z = 2 + °-Z= 3^ ' 



(18) 



(19) 
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These values are the exponents for the GKE equation in 
dimensions d < d c , and are exact to all orders in e due to 
the non-renormalization of A and D in the case r < 2a, 



cf. Eqs. (|12D and (|13|). 

Next, we discuss the case t > 2a. We again remark 
that for any r > 2a we get the behavior for t = 2a, and 
therefore we only have to discuss the latter case |l7j . The 
critical dimension is d c = 2 — a. With e = d c — d, the 
nontrivial FP becomes for a ^ i 



9 



2(2e + 2cr-l) 2(2(7-1) 



0(e 2 



(20) 



The case < a < i; Here, g* is negative, and a FP ex- 
pansion in powers of e does not exist. The KPZ equation 
is a well-known example of this failure of the e expansion. 
The two-loop results for the KPZ equation also show the 
failure of the e expansion around the critical dimension 
4 = 2@. In Fig. | we show the RG flow for g in the 
case where a < \. 

In order to obtain the exponents we can use the one- 
loop result ( |l7| ) of the g* fixed point. For the KPZ equa- 
tion this gives the exact exponents in d = 1 (cf. [Q), but 
despite this "success" the method is uncontrolled due to 
the fact that g* , or A*, is not small at the FP, which 
has been the underlying assumption under the whole RG 
calculation and series expansion in A. The direct sub- 
stitution of ([It]) into the expressions for the exponents 
results in the values (with r — 2a) 



and 



(2 - d){2 -a -d) 
2(3 - 2d) ' 



(2-d)(2-a -d) 



2(3 - 2d) 



(21) 



(22) 
a)/2. 



which in d = 1 yields a = (1 — a)/2 and z = 3(1 
For a = this reduces to the KPZ exponents M. 

The case a > \: In Fig. ^ we show the RG flow for 
g in this case. For dimensions d > | > d c the g* fixed 
point is positive; probably this is an artifact due to the 
one-loop result. Now the e expansion is possible. We can 
obtain the exponents to first order in e at the 0(e) fixed 
point ( po|) with the result 



and it corresponds to a phase transition in the model, 
cf. Fig. ||(c). The phase transition is between the 
"smooth" phase and the strong-coupling rough phase. 
When a > \ the g* fixed point is negative, i.e., the e' 
expansion is only meaningful for a < \. The exponents 
associated with the roughening phase transition are to 
first order in e' 



a 



2(1 - 2a) 



(26) 



For a — the exponents reduce to the well-known results 
a = and z = 2 for the KPZ equation (see, e.g., JL8|). 

V. CONCLUSIONS 

In this paper we have performed a renormalization 
group analysis of the GKE growth equation which is an 
equation with a generalized conservation law described 
by an integral kernel. The results for the FPs and criti- 
cal exponents have shown that the GKE equation encom- 
passes a range of different universality classes, cf. Fig. ||. 

For r < 2a, every point represents a distinct univer- 
sality class, with the SGG and MBE models belonging to 
this case. For all these universality classes we were able 
to obtain the critical exponents exactly, and the values 
are given m Eqs. (|l|) and_([l9|). Furthermore, the expo- 
nents fulfill the identity (Q). 

For t > 2a, every vertical line represents a different 
universality class. The KPZ equation belongs to this 
case. Moreover, we noted the breakdown of the e ex- 
pansion for a < i. As a consequence, estimates for the 
critical exponents could only be obtained by a direct sub- 
stitution of the g* FP value into the expressions for the 
exponents, resulting in the values in Eqs. ( pl| ) and (|2^). 
However, for a > h the e expansion could be used to ob- 
tain the exponent values as given in Eqs. ( p3| ) and (|2~i|). 

In the GKE equation we have the two free parameters 
r and a. As a result, we can always choose these val- 
ues in order to obtain agreement with values for a and 
z determined in an experiment. However, unless one can 
argue that the experiment contains non-local growth ef- 
fects described by a kernel this does not give the true 
explanation of the experiment. 



2(2ct - 1) 



and 



2(2cr- 1) 



0(e 2 



e + 0(e 2 ). 



(23) 



(24) 



For t = 2(7 (t > 2er) there is the possibility of carrying 
out an expansion in the quantity e' = d — d c = d + a — 2 
for dimensions above d r - The FP is 



9*{e) = 



1 



1 



2(2e' + l-2cr) 2(l-2cr) 



e' + 0(e' 2 ) (25) 
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FIG. 1. (a) Diagrammatic representation of the GKE equa- 
tion (^). (b) The vertex A which includes integration over 
(q, Q). The q ■ (k — q) is associated with the outgoing mo- 
menta; J = / A jf^a-fr- (c) The contracted noise 2Dfc T , 
from Eq. @. 
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\ d c + a < d 











FIG. 2. Coupling constant flow for r < 2a. 
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FIG. 3. Coupling constant flow for r > 2a, and a < |. 
For dimensions d > d c (c) shows the possibility of a phase 
transition at g* (see the text). 
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FIG. 4. Coupling constant flow for r > 2a, and a > i. 



d c =2-G 




FIG. 5. Universality classes and critical dimensions for the 
GKE equation. Below the line d c — 2 — a, i.e., for r < 2a, 
every point represents a distinct universality class. Above the 
line d c = 2 — a, i.e., for r > 2a, every vertical line represents a 
different universality class. The KPZ, SGG and MBE models 
are shown with solid circles. The circle at a = % > T = 1 divides 
the t — 2a line into two parts. The part with a < 5 where 
the e expansion does not exist, and the part a > | where the 
e expansion does exist. 
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